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We present a derivation of the Markovian master equation by the renormalization group 
method. Starting from a naive perturbative solution of the von Neumann equation, the re¬ 
duced density matrix with the coarse grained time steps is obtained using the assumption of 
short correlation time of the bath field. Then by applying the renormalization group method, 
we show that the dependence of the specific initial time on the perturbative solution can be re¬ 
moved and the Markovian semigroup master equation in the Gorini-Kossakowski-Lindblad- 
Sudarshan (GKLS) form is obtained in the weak coupling limit. 


1. Introduction 

The quantum dynamics of an open system cannot in general be represented in terms of 
unitary time evolution. One useful method of investigating the open quantum system is to 
formulate its dynamics by means of an appropriate equation of motion for its density ma¬ 
trix; this is a quantum master equation. As one of the typical examples of the open quantum 
system, the particle detector modeP^ is widely accepted as a tool for exploring the quantum 
nature of the considering system. It is composed of detector variables (system) with inter¬ 
nal states that interact with external quantum fields (bath). This model has been applied to 
phenomena of particle creations in curved spacetimes (Unruh effect and Hawking radiation). 
This model is also used to detect the entanglement of quantum fields by investigating the 
correlation between two independent detectors interacting with the quantum fields .^^HlHSEllII 

The state of the detector is obtained by tracing out the field degrees of freedoms which 
interacts with the detector (reduced density matrix). As a result, the evolution of the system 
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is determined by the dynamieal map &, whieh is not unitary in general. 

unitary evolution , 

pr(0) =p(0)®p5 -> prit) = U{t,Q)[p{Q)®pB\U\t,Q) 



dynamical map 

P(0) -> p{t) = £(Op(0) 

To describe physically allowable processes, the dynamical map should preserve the complete 

positivity and normalization of probabilities in the course of evolution. 

The dynamical map can be represented by the master equation that determines the evolu¬ 
tion of the system. In general, the master equation becomes an integro-differential equation, 
and the state at the specific time depends on the past evolution (non-Markovian nature). How¬ 
ever, it is possible to recover the Markovian property of the master equation by assuming suit¬ 
able time scales: the decay time of the bath correlation function tg is sufficiently shorter than 
the relaxation time scale of the system ts . By combining this assumption of the time scales 
with the so-called secular approximation (rotation wave approximation), which neglects tran¬ 
sitions via system energy nonconserving processes, it can be shown that the resulting master 
equation has the Gorini-Kossakowski-Lindblad-Sudarshan (GKLS) form and preserves the 
trace and complete positivity of the state in the course of evolution.It is known that the 
Markovian master equation, which preserves the trace and complete positivity, should have 
the GKLS form and generates the dynamical semigroup 

As the derivation of a master equation that does not rely on the secular approximation, 
the time coarse graining method was proposed.l^^ElEiniBy introducing a coarse graining time 
scale, which is assumed to be longer than the decay time scale of the bath correlation function, 
it is possible to derive the GKLS master equation without assuming the secular approxima¬ 
tion. The purpose of this paper is to present an alternative derivation of the GKLS master 
equation using the renormalization group method.^^^ESI 

The renormalization group method is applied as a tool of asymptotic analysis of differ¬ 
ential equations .Eni xhe naive perturbative solution of differential equations often yields 
secular terms due to resonance phenomena. The secular terms prevent us from obtaining ap¬ 
proximate but global solutions. The renormalization group method is one of the techniques to 
circumvent the problem. Starting from a naive perturbative expansion, the secular divergence 
is absorbed into constants of integration contained in the zeroth-order solution by the renor¬ 
malization procedure. The renormalized constants obey the renormalization group equation. 
We shortly review this method using the following example. Let us consider the van der Pol 
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equation 


(fx , , dx 


( 1 ) 


( 2 ) 


where 6 is a constant. For a small e, the naive perturbative solution up to the order e is 

x{t) = + c.c H-, 

where A is a constant and nonresonant terms such as those proportional to are not ex¬ 
plicitly written. The perturbation fails beyond a time scale ~ l/e owing to the secular term, 
which grows as oc t. To improve this naive perturbative solution, we introduce an arbitrary 
renormalization point r by splitting t to t - t + r and the renormalized constant A(t) by 
A = A(t) -I- SA(t), SA{0) = 0, where the counter term 6A is chosen to absorb r. Thus, 


dA(T) = A(0) - A(t) = -e^A(T)(l - |A(t)P), (3) 

and the naive perturbative solution can be written as 

x{t) = A{T)e“ -I- 6 ^ ^ ^ A(t)(1 - \A{T)^)e'‘ -l- c.c -l- ■ ■ ■ . (4) 

As the renormalization point t is arbitrary, by equating t and t, we obtain the renormalized 
solution x{t) = A{t)e^' with A{t) being the solution of the following renormalization group 
equation: 

dA 6 o 

— = -A(l - |Ap). (5) 

dt 2 

This is an amplitude modulation equation that represents the slow dynamics of the system. 
The renormalized solution gives an approximated but global solution to Eq. ([T]) up to 0 { 6 ). 

In this paper, starting from a naive perturbative solution of the von Neumann equation, 
we rederive the master equation in the GKLS form with the time coarse graining by the 
renormalization group method. We will show that the time coarse graining by the renormal¬ 
ization procedure naturally defines the semigroup structure for the renormalized constant in 
discretized time steps of the evolution of the reduced density matrix. We use the unit in which 
fi= \ throughout the paper. 


2. Derivation of Markovian master equation 

2.1 Naive solution of the von Neumann equation 

Concerning the derivation of a naive perturbative solution of the von Neumann equation, 
we basically follow the presentation by Benatti et a/..I^The total system is composed of de¬ 
tector variables (system) interacting with quantum scalar fields (bath). The total Hamiltonian 




J. Phys. Soc. Jpn. 


FULL PAPERS 


IS 


H = Hi+H^ + AV, 


( 6 ) 


where Hi is the system Hamiltonian, H^ is the bath Hamiltonian and V is the interaetion 
Hamiltonian with V = where cr^ is the system operator and Oa is the bath field. 

The strength of the interaetion between the system and the bath is determined by the eoupling 
constant A, which is assumed to be small (weak coupling). The total density operator prit) in 
the Schrddinger picture obeys the von Neumann equation 


d 

—pr - -i[Ho + AV,pt], Hq 
at 


Hi + Hi 


We aim to obtain an equation for the reduced density operator for the detector system 

p{t) = TrsIprCO}. 

First, we introduce the interaction picture of the density operator defined by 


(7) 


( 8 ) 


PrCO = Ulit, to)PTit)Uoit, to), Uoit, to) = e 
where to is an arbitrary initial time. Then, pj obeys 

j^pt = -umrPT), m = u^t, to)vuo(t, to). 


— (t-to) 


(9) 


where V{t) is the interaction representation of V. We integrate this equation from to to t. After 
two iterations, the formal solution is given by 


Pr(0 = pTido) - iA 


A- 


iA I dti 


[m),pr(to)] 


fA,. r 


dt 2 [V(ti),[V(t 2 ),pT(to)]] + 0(Al. 


We assume that during the evolution, the state of the total system is factorized 

Pt(0 »p(t)<8>pB, 


( 10 ) 


( 11 ) 


where ps = PbC^o)- This assumption is justified because the interaction between the system 
and the bath is weak and the correlation between the system and the bath can be neglected 
when the bath time scale ts is shorter than the system time scale ts . As the bath contains an 
infinite number of degrees of freedom, the back action of the system on the bath is negligible 
and the evolved total state can be expressed as the product state (11). By taking the trace of 
the perturbative solution (fTOl) with respect to the bath degrees of freedom, the reduced density 
operator for the system is 


pit) = Pik) 


-A^ f dh r'^?2TrB|[y(ti),[y(t2),p(^o)®Pfi(^o)]]), (12) 

Jtn Jt(\ 
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where we have assumed 

TrB{mp{to)®PB{to)] = 0. (13) 

If we do not eonsider the baek aetion of the system on the bath and the bath variable evolves 
freely, the bath operators in the interaction picture are the same as those in the Heisenberg 

picture 0^(0 = Thus, the condition (fT3l) is equivalent to (Oa) = 0 because 

TVB{V{t)p(to) ^pB(to)} = '^d-A(t){^A(t)) = 0, 

A 

where the expectation value of the bath field 0(t) (Heisenberg operator) is evaluated as 

(0(0) = Trs{0(OpB(to)}. 

We rewrite the integral in Eq. (fT^ . Let us define 

Ai2 = Tvb |[y(ti), [Vit2),pito) ®PB(to)]]), (14) 

and decompose the integral as 

r dti r dt2Ai2 = - f dti r dt2{Ai2+A2i) + — f dt\ f dt2(Ai2 - A 21 ). (15) 

Jlo Jto ^ Jlo Jto ^ Jto Jto 

From now on, we omit the symbol (8) for simplicity. The first integral is 


r 

Jtf) Jtn 


dt2TrBi^{mm2),p(fo)PB(to)} - nh)p{k)PB{kmt2) 


( 16 ) 


and the second integral is 


^ f dp f dt2TrB\[V{ti),V{t2)\,p{tQ)pB{tQ)\ 

^ Jto Jto 

= ^ r r dt2Sgn{ti - t2)[TrB(pB(to)V(ti)V(t2)),p{to)Y 
^ J/n 


(17) 


Using V = Xa and introducing the correlation function of the bath field b)^!^ 

GaiA2(6 ~ h) = TrB{OA((tl)OA2(t2)PB(^)} = (d>A|(E)d>A2(l2))5 (18) 

we can rewrite the integral in the right-hand side of Eq. (fT^ as 

p(t)-p(to) = At x{-i[H12, Po]+-ClPo]), At = t-to, po=p(to), (19) 


where we have defined 


H 


12 




dsi 


J nA.t 

dS2 Sgn(5i - S2 )Ga^^{Si - S 2 ) (Ta^Si + to)(TAlS2 + to), (20) 

0 


^2 nAt i^At . 

JJPq\ = TT y, I dS2GA^A2(Sl - S2)(o-A 

ATI2 ^0 ^0 ^ 


(Si + to)poO'Ai(Si + to) 
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- 2 + tQ)o-A^{S 2 + ^o),Po})- ( 21 ) 

In Eqs. (I20l) and (l2T]) . the system variable (Ta( 0 in the interaetion pieture is replaeed by (TaCO 
in the Heisenberg picture whose evolutions are determined by the free system Hamiltonian 
//q. This is justified within perturbation up to O(A^). The time dependence of the system 
variable is determined by the free Hamiltonian as 

ctaCO = ^ UABit - to) (Tb, mab(O) = 1 , ( 22 ) 

B 

where the specific form of the function mab( 0 depends on the system Hamiltonian. Using this 
relation, we obtain 


Hn = Z 

B[,B 2 

-^^[Po] = ^ ^51^2 (^B2P0 ^Bi “ 2 > 


Bi,B2 


with 


^2 i^At 

^BiB2 = I dsi I Sgn(5i - 52 )GaiA 

ri„ Jo Jo 


(23) 

(24) 

(Si - S2 ) UaiBi(Si)Ua 2 B 2 (^ 2 ), (25) 


c 


BIB2 


Ai,A2 

i2 r'At 

=Tv r 




f 


dsi dS2GA^2(^l - S2 )UaiB,(Si)Ua2B2(S2)- 


(26) 


As Eq. (fT9l) determines p(t) from the specific initial state at to, we cannot obtain the master 
equation even though if we take the limit At —> 0. To convert (fT9l) to the Markovian master 
equation, we introduce the time coarse graining and assume the short correlation time of the 
bath. Then we apply the renormalization group methocP^Hl to recover the Markovian master 
equation. 


2.2 Time coarse graining and renormalization group 

The naive perturbative solution (fT^ has the following structure: 


pit) - pito) 

2 _^ p^t p, 

= A^At X — I dsi I 

am,,B2'^o Jo 


= A^Atxnipo], 


dS2GAMSi - S2 )UaiBi(Si)Ua2B2(^2)^BiB2[^1 “ ^2,Po] 

(27) 


where FbiB2 represents a combination of the following terms: 

Sgn( 5 'i - S2)crBi(TB2P0, Sgn( 5 i - 5 ' 2 )PoC’'B|CrB 2 , PoO"B|C’'B2> (^B2PoO'Bi- 
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Let us introduce a coarse graining time scale A and consider the evolution of the state by this 
time step. Using (1271) . the evolution of n time steps from to to to + ^A is 


p(to + nA) - p(to) 

Y ^nA ^nA 

= A^A X — I dsi I 

A Jo Jo 

= d^Ax2)„[po]. 


dso 


^ - S2)UAiBi(Sl)UA2B2iS2)FBiB2[Sl “ 52,Po] 


Ai,A2,Bi,B2 


(28) 


At this point, we assume that the coarse graining step size A is sufficiently larger than the cor¬ 
relation time Ib of the bath field and the relation tg A holds. By introducing new integration 
variables as ti = si - S 2 ,T 2 = 52 - ■s'l, the integral of (1281) is 


X niA ^n/A 

dTi J dS2 E G{Ti)u{Ti + S2)u{S2)F[Ti] 

X nA pnA 

dT 2 J ds\ E G{-T2)u{Si)u{T2 + 5l)F[-r2] 

X oo ^nA 

dr J ds^(G(T)u(s)u(s)F[T] + G(-t)u(s)u(s)F[-t]) , 


(29) 


where we have extended the upper bound of the Ti,T 2 integrals to infinity because the corre¬ 
lation function decays rapidly for ti , T 2 ^ 0 and the main contribution to the integrals comes 
from the region ti,T 2 » 0. Then, !D„ is reduced to the following form: 


1 V—1 F^ 

F>„[po] = - ds > UA^Bi{s)UA 2 B 2 iA) 


dTGAiA2(T)FBiB2[FPo], 


(30) 


and it is possible to derive the following relations for n > 1: 

^ ^nA _^ ^+oo 

(D„ - 2)„_i)[po] = - I dsy^ u(s)u(s) I dr G(t)F[t, po] 

^ J(n-1)A J-oo 

-A 


-U'-S 


u{s + (n - 1)A) u(s + (n - 1)A) 


f 


dr G(t)F[t,po] 




[ 


^-iHl («-l)Ap^ ^iHl («-l)Aj (n-l)A^ 


(31) 


where we have defined Do = 0. Let us denote the state at to + nA as Pn = p(to + nA). 

The perturbative solution (ITS!) is represented as 

p„ = C + (d2A)D„[C], C=po=p(to). (32) 

As the state in the interaction picture does not evolve in 0{A^), C is a constant of motion in 
this order. Hence, we regard the 0{A^) term in (l32l) as the secular term due to the interaction 
and apply the renormalization group method to remove the secular behaviour of the solution. 
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We introduce an arbitrary renormalization point fc(0 < k < n). The renormalized constant Ck 
and the counter term 6Ck are introduced as C = Q + 6Ck, Co = C, 6Co = 0. The counter term 
is chosen so as to absorb the secular term of the solution: 

6Ck = Co-Ct = -(A^A)!DdCo]. (33) 

By using the renormalized constant, the naive perturbative solution (l3^ can be written as 

p„ = Ck + (d'A)(©„ - D,)[Co], (0<k< n). (34) 

Equation (13^ defines a map : Cq i—> Ck (0 < k). As the renormalization point is arbitrary, 
the following two equations hold for two different renormalization points ki , ^2 (0 < ki,k2 < 
n): 

Ck, - Co = {A^^){!Dk, - Do)[Co], C,, - Co = {A^^){!Dk, - ©o)[Co]. 

By subtracting these equations side by side, we obtain the following relation up to O(A^): 

Ck, - Ck, = (A^A) (Dk, - Dk,) [Ck,], (0<k:< ^ 2 ). (35) 

Thus, this equation defines a map 'Rk,-k, '■ Ck, 1 —> Ck, (0 < k 2 - ki). As the map defined 
by (l35l) satisfies the composition law 'Rk,-k, ° ^k, = 'f^k, ih < ^ 2 ), it generates a discretized 
version of the dynamical semigroup. Therefore, the evolution of the system can be Markovian 
in the discretized (coarse grained) time steps. In Eq. (I34l) . as the renormalization point k is 
arbitrary, we can choose k = n: 

Pn = Cn. (36) 

This is the renormalized solution and the renormalized constant C„ satisfies the discretized 
version of the renormalization group equation (l35l) . Combining (l35l) with (l36l) and using the 
relation (l3T]) . we obtain the following difference equation for p„: 

pn-pn-l = (d^A)(£)„ - 2)„_i)[p„_i] 

= (T^A) (37) 

Now, we introduce a new time variable r = A^t and take the weak coupling limit T —> 0. Eor 
p(t), we define the time coarse grained state as 

p^^iTn) = P(t„//12)|^^q = Pnl^^o . (38) 

To determine the evolution for the finite interval of time t„-to = A^nA, we must take « —> oo. 
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In this limit, 


and 


lim 


Pn Pn— 1 


,. P^^iTn) - P^^{Tn “ A) 

- 



(39) 


limp„_i = limp‘^°(T„ - A^A) = p^^(t„), (n - 1)A = lim \-{t„ - d^A - tq) = (t„ - Tq)IA^. 

/l-»0 /i^o A^O A^ 

Therefore, we obtain the differential equation for the state p^^{t)\ 

dr ( J 

The differenee equation (ITTI) is eonverted to the differential equation in the resealed time 
variable r with the small eoupling limit. Henee, (|4^ represents the evolution of the time 
eoarse grained state. 

Returning baek to the Sehrodinger representation using p^’^it) = 
we obtain the following linear Markovian master equation: 

^pC^(0 = -i[Hlp^%)] + D,[p^%)], (41) 

where the bath-dependent eontribution Dj \sf^{t)\ eontains both a Hamiltonian and a dissi¬ 
pative term 




(42) 


with 

^ BuBi 

-C^[p] = Yj ^BiB2 (^B2P (^Bi - ^ {crB,(rB2,P^ ■ (44) 

B,,S 2 ^ ’ 

Thus, the master equation is 

= -ml + Hl^,p^^{t)] + £^''(0]. (45) 

The time-independent eoeffieients ^BiB 2 given by (1251) and (l2^ replaeing At —> A. 
These eoeffieients do not eontain the initial time and the initial state and depend only on the 
eoarse graining time A. Therefore, if the eoeffieients Cgjg^ ^ positive matrix, the master 
equation (l4T]) has the GKLS form and the state evolves preserving the traee and eomplete 
positivity. To show the positivity of the matrix defined by (26), we must eheek the 
inequality £ 5 , 62 Tb, ^b^B 2 Xb 2 > 0 for an arbitrary veetor;\fB. The form of the matrix is 

^b,B2 - ^ dSi(^Ai{Si)UAiBiiSi) u dS20A2iS2)UA2B2iS2)^ 


9m 
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_A\ , 

where we introdueed the operator 




■?r 




(46) 


(47) 


Without loss of generality, we ean assume that the operators Oa and cta are Hermitian. Then 
the funetion UAsit - h) representing the evolution of the system variable cta is real and the 
operator La is the Hermitian Lb = L\. Hence, 

A 




B1B2 


42 


B\,B2 


A 

> 0 , 


^T / 




V Bi 


' \ B2 


(48) 


because this quantity is an expectation value of the positive operator. Therefore, the coef¬ 
ficients form a positive matrix. The master equation (l45l) is the same one derived in 
Refs. 5 and 12-14. 


3. Summary 

In this paper, by the renormalization group method, we have rederived the Markovian 
master equation in the GKLS form under the following assumptions: 

(1) The interaction between the system and the bath is weak (weak coupling) and the back 
action of the system on the bath is negligible. 

(2) The correlation time of the bath field is sufficiently shorter than the relaxiation time of 
the system and the existence of the factorized state at the specific initial time to can be 
assumed. 

(3) The coarse graining time A is longer than the time scale of the bath field. 

The naive perturbative solution Eq. (1271) does not have a form of the master equation. To 
transform it to the master equation, we introduced the coase graining time scale A, which is 
larger than the bath time scale tg. This reduces the double integral with respect to time in 
the solution to the single integral, and this reduction is equivalent to introducing the Marko¬ 
vian approximation carried out in conventional derivations of the master equation. However, 
at this stage, the equation does not have the form of the master equation. Then, we applied 
the renormalization group method and eliminated a specific initial time. Conventional deriva¬ 
tions of the master equation usually do the same thing by requiring the time translational 


10fl2] 
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invariance of the eonsidering system and the strueture of the equation does not depend on 
the speeific initial time. We have carried out the equivalent procedure by assuming that the 
strueture of the evolution does not depend on the speeifie renormalization point and derived 
the renormalization group equation, which provides the desired master equation. 

The obtained renormalization group equation generates a dynamieal semigroup for the 
renormalized eonstant, and the renormalized density operator for the slow time variable r 
satisfies the time loeal Markovian master equation in Lindblad form. As diseussed in Refs. 5 
and 12-14, the master equation (l45l) reduces to the one with the rotational wave approximation 
and the seeular approximation in the limit A ^ oo. However, for the purpose of deteetion of 
the entanglement of the quantum field using particle deteetors, the master equation with the 
time eoarse graining is more suitable beeause it ean deteet the vaeuum quantum fluetuation, 
which is the cause of the entanglement of the quantum field. 
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